First Experimental Observation of Superscars in a Pseudointegrable Barrier Billiard 
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With a perturbation body technique intensity distributions of the electric field strength in a flat 
microwave billiard with a barrier inside up to mode numbers as large as about 700 were measured. 
A method for the reconstruction of the amplitudes and phases of the electric field strength from 
those intensity distributions has been developed. Recently predicted superscars have been identified 
experimentally and — using the well known analogy between the electric field strength and the 
quantum mechanical wave function in a two-dimensional microwave billiard — their properties 
determined. 



PACS numbers: 05.45.Mt, 03.65.Sq, 42.25.Fx 

Planar polygonal billiards with angles Oij — irnij/nj, 
where rrij and rij are co-prime integers, have been stud- 
ied both classically and quantum mechanically 1] . When 
7ti j =/= 1 the motion in phase space is not restricted to a 
torus like for integrable systems, but to a surface with 
a more complicated topology. Accordingly, such planar 
polygonal billiards are called pseudointegrable. It was 
established numerically (see @, Q and refs. therein) that 
the statistical properties of the eigenvalues of the cor- 
responding quantum systems are intermediate between 
those of a regular and a chaotic system. The properties 
of the wave functions of planar polygonal billiards are 
also intriguing, as they show a strong scarring behavior 
which can be related to families of periodic orbits 0] • In 
a plot of eigenfunctions in the barrier billiard scars are 
clearly distinguishable from non-scarred eigenfunctions. 
This pronounced scar structure does not disappear at 
large quantum numbers in contrast to that in chaotic 
systems 0, H, 0| • To stress this difference it is proposed 
in 0] to call the scars in pseudointegrable systems super- 
scars, an expression used by Heller in his early seminal 
paper in a different context. The aim of this letter is 
to report on the experimental investigation of superscars 
in the barrier billiard. This is a rectangular billiard of 
area l x ■ l y which contains an infinitely thin barrier. In 
the experiment presented here the barrier is placed on 
the symmetry line x — l x /2 and its length equals l y /2, 
where l y is the length of the shorter side of the rectangle. 

In the present work the quantum barrier billiard is sim- 
ulated by means of a microwave billiard. Microwave bil- 
liards are flat cylindrical resonators @, 0] . Below the crit- 
ical frequency f c = c/2h, where c is the velocity of light 
and h is the height of the cavity, the electric field is the 
solution of the scalar Helmholtz equation with Dirichlct 
boundary conditions. This equation is mathematically 
equivalent to the Schrddinger equation for a quantum 
billiard of corresponding shape (see e.g. @, @|). 

The experimental observation of the superscars pre- 
dicted in [4| requires (i) the detection of eigenmodes at 



sufficiently high quantum numbers N and (ii) the reso- 
lution of all modes up to these quantum numbers. The 
number N(f c ) of resonances below the critical frequency 
f c is approximately given by the first term of the Weyl 
formula 0, i.e. N(f c ) = & (2p / c ) 2 - £ (f f , where 
A is the area of the billiard. Hence, the number of exper- 
imentally accessible resonances increases with the area 
and decreases with the height of the billiard. The reso- 
nances can be well resolved, if their widths are small in 
comparison with the average spacing between adjacent 
resonances, i.e. a microwave cavity with a high quality 
factor Q is needed. The latter is proportional to the ra- 
tio of the volume to the surface of the cavity, therefore for 
flat cylindrical resonators to the height of the resonator. 
A compromise between a high Q and a large number of 
resonances is obtained by designing a rectangular cavity 
with side lengths of l x = 1000 mm and l y — 585 mm, and 
a height of h = 25 mm. The barrier has length 297 mm, 
a width of 4 mm and the same height as the rectangular 
cavity. As in our earlier experiments (llj the resonator 
consists of a bottom, a lid and a frame in between. All 
parts are squeezed together with screws. The bottom 
and lid plates respectively are manufactured from pol- 
ished copper. The frame has the shape of a rectangle and 
consists of 10 brass segments. To ensure a proper electri- 
cal contact even at high frequencies, wires of solder are 
placed between the bottom, the frame and the lid. The 
so constructed microwave barrier billiard has about 700 
resonances below f c = 6 GHz, the average quality factor 
is Q = 1.8 • 10 4 and the mean resonance width is at least 
10 times smaller than the mean resonance spacing. 

In order to obtain information on the wave functions 
in the quantum barrier billiard we determined the elec- 
tric field strength distribution in the corresponding mi- 
crowave billiard. Thus first intensity distributions of the 
electric field strength were measured with the so called 
perturbation body method [13]. A small metallic body 
alters the resonance frequency fo of the cavity, where ac- 
cording to the Maier-Slater theorem [l3| the frequency 
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shift depends on the difference of the squares of the elec- 
tric field strength E and the magnetic field B at its lo- 
cation inside the cavity. 



A/ = f ( Cl E 2 -c 2 B 2 ). 



(1) 



Here, c\ and ci are constants determined by the geom- 
etry and material of the perturbation body. In order to 
obtain the electric field strength distribution, the per- 
turber needs to be moved inside the closed resonator. 
This is best done with an external guiding magnet. How- 
ever, since we are interested in the electric field strength 
only, the magnetic field contribution in Eq . (TO has to be 
removed. In earlier experiments [12, EH, EH EH with 
microwave billiards the magnetic field contribution to 
the frequency shift was minimized by using a needle like 
metallic bead. We, for the first time, used a nonmetal- 
lic perturber (of so called magnetic rubber) which is a 
rubber- like plastic combined with barium ferrite powder, 
whose grain size is of the same order of magnitude (1-10 
/Ltm) as the skin depth of radio frequency waves. There- 
fore eddy currents cannot propagate, such that the mag- 
netic rubber does not interact with the microwave mag- 
netic field B. Accordingly there is no contribution to A/ 
in Eq. ©. 

We first performed test measurements with the rect- 
angular cavity and a cylindrical perturbation body with 
a height of 5 mm and a diameter of 2 mm. Typically 
the positive frequency shift induced by the electric field 
strength is of the order of 100 kHz. As expected, we did 
not detect any negative frequency shifts, which would re- 
sult from the magnetic field term in Eq. |T]) within the 
experimental uncertainty of the frequency shift measure- 
ment of about 1 kHz. Moving the perturbation body 
with an external guiding magnet across the billiard sur- 
face with a spatial resolution of about one-tenth of a 
wavelength and measuring the frequency shifts A/ at 
each of its positions yields the intensity distribution of 
the electric field strength in the whole microwave bil- 
liard. The shift Af is determined as in [17]] via phase 
shift measurements. In the frequency list mode the vecto- 
rial network analyzer HP8510C allowed the simultaneous 
measurement of the phases at up to 30 different frequen- 
cies. This enabled us to crucially reduce the measure- 
ment time, we measure 30 intensity distributions in one 
raster scan for a grid with 200x100 points in 12 hours. 
To minimize the influence of a temperature drift, which 
causes an additional frequency shift of about 40 kHz/K, 
the temperature of the microwave billiard was stabilized 
to within 0.1 K. 

Due to the symmetry of the barrier billiard there are 
two types of modes, anti-symmetric and symmetric ones 
with respect to reflections at the barrier, respectively. 
The former are trivial as they coincide with those of the 
rectangular billiard, the latter are exactly those we are 
interested in. We were able to measure, with the ad- 
vanced experimental method described above, the inten- 



sity distributions of altogether 590 modes for level num- 
bers N =90-680. In contrast to the technique described 
in [HI with the perturbation body method we detect the 
E field intensity only. For a detailed investigation of 
the symmetric mode properties we developed a special 
procedure to reconstruct the E field, i.e. the quantum 
mechanical wave function tp, itself and its sign from the 
measured intensity distribution. 

For this we fit the square of a sine series to the mea- 
sured field intensity. Along a line of constant x, x = xq, 
in the billiard plane the latter can be represented as 



tpl{x = x ,y) = 



-i 2 



]Ai(x ) sin 



(2) 



where ki — wi/ly, l y is the length of the billiard along 
the line x — xq, and Ai{xo) are the expansion coeffi- 
cients. With this choice of ki the Dirichlet boundary 
condition is automatically fulfilled. The number of terms 
giving a significant contribution to the sum can be esti- 
mated semi-classically, i.e. n ~ 2l y /X and A = 2n/k. For 
the determination of the expansion coefficients Ai(xo) 
in Eq. ([2]) a nonlinear fit routine which is based on the 
Levenberg-Marquardt algorithm [l9| was used. The non- 
linear fit requires a trial vector of coefficients as initial 
values. We used random numbers which are uniformly 
distributed in the interval from -max(| , 0(x = xo,y)\) to 
max{\^>(x = xo,y)\) as entries of the trial vector, whose 
typical dimension is about n — 20. The nonlinear fit 
has no unique solution as the Levenberg-Marquardt al- 
gorithm only provides a set of coefficients Ai(xo) which 
gives a local minimum in \ 2 ■ To find the global mini- 
mum we generated about 1000 trial vectors and carried 
out a fit for each of these. Then that set of coefficients 
Ai(x$), which gives the smallest x 2 > provides the wave 
function values along the line with constant x — xq . Due 
to the continuity of the wave function for the neighboring 
lines we can use the so determined coefficients as entries 
of their trial vectors. An example for the measured in- 
tensity and the reconstructed field of the eigenmode with 
N = 613 and frequency / = 5.48004 GHz is shown, re- 
spectively, in Fig.[TJ From 290 measured intensity distri- 




FIG. 1: (Color online) (a) The measured intensity distribution 
for the mode N — 613. (b) The reconstructed E field, i.e. the 
wave function. The corresponding color scales are at r.h.s. of 
the graphs. 

butions corresponding to the symmetric wave functions 
we reconstructed 230. For the remaining 60 distributions, 
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due to large noise or nearly overlapping resonances, re- 
construction was not possible. 

Among the large number of wave functions several su- 
perscarring ones were observed. Four examples of de- 
tected superscars are shown in Figs. [2] and [3l Each wave 




FIG. 2: (Color online) The top row shows examples of ex- 
perimentally obtained superscars in the barrier billiard. The 
color scale is the same as in Fig. [lib). l n the bottom row the 
corresponding classical orbits are indicated (dashed lines). 




FIG. 3: (Color online) The same as in Fig.[2]but for two other 
types of superscars. 

function shows a clear structure which is obviously re- 
lated to the corresponding periodic orbit (PO). 

To understand the quantum nature of the superscars 
we briefly outline how the PO families of a barrier bil- 
liard can be constructed. For this purpose we mark one 
corner of the barrier billiard with a black point as shown 
in Fig. 0Ja), thus defining the initial orientation of the 
billiard. We follow the orbit of the propagating particle 
(dashed line in Fig. IUa)) and mirror the billiard each 
time the particle hits a side of the billiard. This proce- 
dure is repeated until the mirrored copy of the billiard 
resumes its initial orientation. All parallel trajectories 
with the same length form a family of POs. Each family 
corresponds to an infinitely long periodic orbit channel 
(POC). This POC is confined to the region between the 



two straight lines which connect the singularities, i.e. the 
barrier tip and its images. These lines are called the sin- 
gular diagonals (SD). In the semi-classical limit we can 
associate with each classical trajectory a wave. Waves 




FIG. 4: (Color online) (a) Illustration of the construction of 
the PO channel. The black point in the upper left corner 
of the original billiard fixes its initial orientation, (b) The 
unfolding of the experimentally obtained superscarring wave 
function with number N = 587. It has 110 wave maxima 
along the POC and one perpendicular to it, corresponding to 
the quantum numbers m = 55 and n = 1, respectively. 

travelling in a POC are scattered at an infinite array of 
images of the barrier tip. At the SD the intensity of the 
scattered wave is proportional to where k is the 

wave momentum [20|. Therefore, it tends to zero and 
fulfills the Dirichlet boundary condition there with in- 
creasing k. Accordingly, in the semi-classical limit the 
SD act as perfect mirrors such that in each POC quasi 
states can be constructed. These are called unfolded scar 
states in Q and are given as 

= z~ lkmi ■ «n (—V) ■ X(V)- (3) 
V w / 

Here £ is the coordinate along the POC and ij (0 < r\ < 
w) the one perpendicular to it, and w is the channel 
width. The quantity xiv) is t ne characteristic function 
of the POC (x(rf) = 1 if < 77 < w, and x(v) = 
otherwise). The frequency of the eigenfunction (Eq. pj)) 
of the POC is 

'"•" = si/* + (t)'- (4) 

Due to the Dirichlet boundary condition the wave gains 
a phase of tt each time the classical trajectory crosses the 
billiard boundary. If the PO experiences an even number 
of boundary crossings while travelling through the POC, 
then k m = 2nm/l, otherwise k m — tt (2m + l)/l, were I 
is the length of the PO. When folding back the unfolded 
superscar states of Eq. ([3]), complicated expressions for 
folded scar states >5 ' m .n(x,y) are obtained which, except 
on the SD, obey the Schrodinger equation for the barrier 

billiard (A+ (^/ m ,„) 2 ) * m , n (x,y) = with Dirichlet 
boundary conditions on the billiard boundary. As can 
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be seen in Fig. [4|b) the experimentally obtained super- 
scarring wave functions are mainly concentrated inside 
the corresponding POC. This confirms the approximate 
validity of the ansatz for the superscar wave functions 
given in Eq. ([3]). 

A typical eigenfunction of the barrier billiard may have 
contributions from many scar states. In order to ob- 
tain a quantitative measure for this we followed Ref. Q 
and computed the overlap integral of folded scar states 
y) given in Eq. ([3]) with the measured wave func- 
tion tyz(x,y) at the unfolded resonance frequencies fx, 
where the unfolding was performed as usual 8, 9] using 
the Weyl formula 0, 

C m , n (h) = J*%Z(x,y)* h (x,y)6xdy. (5) 

In the computations we fixed the number of excitations 
perpendicular to the POC to n = 1; m is obtained from 
the condition that \fx — fm,n\ should be minimal. The re- 
sulting squares of overlap coefficients |C mjn (/)| 2 , versus 
the unfolded resonance frequency / for two different fam- 
ilies of superscars, are shown in Fig.O The frequencies of 
the superscars are given by Eq. {4|. Near almost all pre- 
dicted superscar frequencies, marked by points in Fig. [5j 
we found a wave function with a strong contribution to 
the integral of overlap with a constructed superscar state. 



different m, 

Pn(f) = (j2\C m Ah)\ 2 S(f-h + fm,n)) ■ (6) 
\ x ' m 

In the bottom row of Fig. [5] this quantity is shown for the 
two scar families in the top row of the figure. It is evident 
that the constructed superscars Eq. © have a certain 
overlap with the neighboring states. The finite width 
of the local density curves indicates that the Dirichlct 
boundary condition along the SD is an approximation, 
that is, the field can leak out of the POC. As can be seen 
in Fig. [5j it seems to be reasonable to approximate the 
averaged local density of the superscar states by a Breit- 
Wigner shape, but further work on this is in progress. 

In summary, we have measured the electric field 
strength distributions at 590 resonance frequencies of 
the barrier billiard and reconstructed almost all corre- 
sponding symmetric wave functions. Many of the wave 
functions show a surprisingly pronounced structure of 
so called superscars which is connected with classical 
PO families. We plan to extend our measurements of 
wave functions in barrier billiards with one or more non- 
symmetrically placed barriers. The statistical properties 
of the nodal domains of polygonal billiards are another 
very interesting problem, which will be investigated. 

This work has been supported within the DFG grant 
SFB634. 
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FIG. 5: The top row shows the overlap of measured wave 
functions of the barrier billiard with analytically constructed 
scar states for different POs. The points indicate the pre- 
dicted positions of the superscar at the rescaled frequency /. 
In the bottom row the corresponding local density is plotted. 
The solid lines indicate the best fit of a Breit-Wigner shape 
to the data. The insets show the corresponding PO families. 



To analyze the overlap of the constructed superscars 
with the neighboring non-scarred states, we calculated 
the average over all the spikes showed in Fig. and ob- 
tained the so-called local density of states averaged over 
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